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Abstract 

In [H E] it has been shown that galactic potentials can be kinematically linked to the 
observed red/blue shifts of the corresponding galactic rotation curves under a minimal set 
of assumptions: the emitted photons come from stable timelike circular geodesic orbits of 
stars in a static spherically symmetric gravitational field, and propagate to us along null 
geodesies. It is remarkable that this relation can be established without appealing at all to a 
concrete theory of gravitational interaction. Here we generalize this kinematical spherically 
symmetric approach to the galactic rotation curves problem to the stationary axisymmetric 
realm since this is precisely the symmetry that spiral galaxies possess. Thus, by making use 
of the most general stationary axisymmetric metric, we also consider stable circular orbits 
of stars that emit signals which travel to a distant observer along null geodesies and ex- 
press the galactic red/blue shifts in terms of three arbitrary metric functions, clarifying the 
contribution of the rotation as well as the dragging of the gravitational field. This station- 
ary axisymmetric approach distinguishes between red and blue shifts emitted by circularly 
orbiting receding and approaching stars, respectively, even when they are considered with 
respect to the center of a spiral galaxy, indicating the need of precise measurements in order 
to confront predictions with observations. We also point out the difficulties one encounters 
in the attempt of determining the metric functions from observations and list some possible 
strategies to overcome them. 
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1 Introduction 



The galactic rotation curves provide a direct method of determining the gravitational field 
inside a spiral galaxy since they have been measured for a great amount of galaxies HI [5] . 
These curves are obtained by measuring the red/blue shifts of light emitted from stars and 
from the 21 cm radiation from neutral gas clouds. The observations show evidence that the 
red/blue shifts z, or equivalently, the tangential velocities of rotation v, remain constant or 
decay more slowly than the Keplerian behaviour (v 2 ~ 1/r) up to distances far beyond the 
luminous radius of these galaxies. By performing a naive Newtonian analysis of this effect, 
one deduces that the energy density of the galaxies decreases approximately as r~ 2 , and 
hence, the mass of these bodies should increase as m(r) ~ r. Since the observed luminous 
galactic components do not produce this growing behaviour, a question arises: what is the 
reason of such an effect? Nowadays there is a strong believe that dark matter is responsible 
for it, being the major bounded constituent of galaxies and galaxy clusters (~ 25% of the 
total energy density of our Universe [6]). However, one can naturally ask whether this large 
unseen mass does not produce a relevant gravitational red shift. Alternative approaches like 
modifications of newtonian dynamics based on these observations have also been developed 
[7]-|ll|. Here, following [I], we shall analyze the problem on the basis of what it is directly 
observed: the red/blue shifts. This approach enables us to keep track of the effect of the 
underlying made assumptions, and to be aware of when they are not longer valid. 

There have been previous approaches to the galactic rotation curves problem that make 
use of stationary axisymmetric metrics within the framework of General Relativity (see the 
comoving dust solution of Cooperstock et al. [12j [T3] ) . Within these models it is claimed 
that the known data from galactic rotation curve can be described with at most relatively 
little extra matter by non- linear general relativistic effects in galactic dynamics. An im- 
portant point within their approach is that even though fields and velocities are small in a 
galaxy, it is not consistent to describe the latter through a Newtonian approximation since 
there are non-linear contributions of non- negligible size coming from Einstein equations. 
However, these models have received criticisms in several directions that point out to un- 
physical features like the need of additional exotic matter source in the galactic disk |14| 
or infinite mass at large distances |15j . the presence of singularities when continuing the 
interior solution into a consistent exterior configuration [16j . among others, as well as the 
inconsistency of the use of comoving frames with the condition of differential rotation |17j . 
An important observation upon this model was pointed out in |18j . namely, the stationary 
axisymmetric metric of the comoving dust solution of Cooperstock et al. does not possess 
the most general form since it has only three arbitrary functions instead of four, moreover, 
their assumed Weyl gauge (W = r) is not consistent with the Einstein equations since this 
metric function is not harmonic and hence, it does not belong to the most general class of 
stationary axisymmetric solutions, the Lewis-Papapetrou class, a circumstance that might 
be connected with the problems of the model. 

The aim of this paper is to provide a stationary axisymmetric kinematical description of 
the galactic rotation curves problem since this is a more realistic symmetry compared to the 
spherical one when studying spiral galaxies. We make use of the most general stationary 
axisymmetric metric and express the galactic red/blue shifts measured by a distant observer, 
that can in principle be compared to observations provided by astronomers, in terms of three 
arbitrary metric functions. We also clarify the contribution of the rotation and the dragging 
of inertial frames due to the gravitational field. We further point out the difficulties we have 
when determining the metric functions from observations without making reference to any 
theory of gravitational interactions, and comment on some possible ways to overcome them. 
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2 Stationary Axisymmetric Rotation curves 



We shall start by assuming that stars behave like test particles which follow time-like 
geodesies of a rotating axially symmetric space-time associated with galactic disks which 
possess such a symmetry. The most general line element for a space-time of this kind has 
the following form: 



ds 2 = -e^dt 2 + Q 2 dr 2 + R 2 d9 2 + sin 2 9 (dtp- Wdt) 2 , (1) 

where &,Q,R and W are all functions of r and 9. We shall also consider two observers Oe 
and Od with 4- velocities u E , u D , respectively. Observer Oe corresponds to the light emitter 
(i.e., to the stars placed at a point Pe of space-time), and Od represents the detector at 
point Pjj, which is located far away from the light emitter and is ideally located at r — > oo. 

We further assume that stars move on the galactic plane and, thus, the polar angle can 
be fixed 9 = vr/2, so that u% = (f7*, U r , 0, U^) E , where = x» and the dot stands for 
derivation with respect to the proper time of the particle (star). 

On the other hand, the 4-velocity of the detector located "far away" from the source, in 
a stationary axisymmetric background, is given by u D = (Z7*, U r , 0, U v )d in the language 
of the above mentioned coordinates and conventions. In this expression the component 
of the 4-velocity u D accounts for the dragging of the observer at Pd due to the rotation of 
the galaxy. This effect must be taken into account when considering the measurement of 
red/blue shifts of light signals emitted in our own galaxy or in galaxies "close" to ours, for 
instance. However, if the studied red/blue shifts correspond to galaxies located "far away" 
from us, so that we can neglect the dragging effect, then we can consider that the detector 
is static, i.e., that the Od's 4-velocity is tangent to the static Killing field JL and hence its 
4-velocity is given by u D = (U t , U r , 0, 0)d- Later on we shall give a qualitative estimation 
of whether we can consider this ideal limit (when the dragging effect can be neglected) in 
terms of the contribution of the angular velocity of a given galaxy to the measured red /blue 
shift. 

We further normalize the 4-velocity as usual (u^u^ = —1), a condition which renders 
the following relation 

-l=g tt {(J*) 2 + g rr (U r ) 2 + g w {U*f + ^t^U* . (2) 

The stationary axisymmetric metric ([T]) possesses two commuting Killing vectors: the time- 
like = (1,0,0,0) and the rotational one tp^ = (0,0,0,1). The corresponding conserved 
quantities are the energy and the angular momentum per unit of mass at rest of the test 
particle and read 

E = -g^v? = - (gull* + g ttp m) , (3) 

L = g^u u = g w U^ + g^U* . (4) 

These relations are useful to obtain the expressions for [7* and U v in terms of the metric 
and these conserved quantities: 

TT t _ E 9w + Lgup _ E-LW 
9up ~ 9ttg w e zq > 

_ Eg tv + Lg tt _ L (E - LW) W 
9t<p ~ 9ttg w K e 
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By introducing these 4-velocity components ([5]) and © in the line element ([2]) we get 

9rr (U r f = EU l - LU^ - 1 , (7) 
or, after multiplying by e 2 * we equivalently have 



r2 2* 

q2 {ur) 2 + e 2$ + e_ _ {E _ LW)2 = ^ q2 {Jjr)2 + V ^ {E ^ L ^ = Q _ (8) 



This equation resembles an energy conservation law for a non-relativistic particle with 
position dependent mass moving in an effective potential that depends on E and L. This 
equation can be further reduced by considering that the orbits of stars are circular U r = 0: 

t2 2$ 

V ef f = e 2 * + -{E- LWf = . (9) 

On the other hand, the time-like geodesies that follow the stars orbiting around the 
galaxies can be written as 

dTJV 

— + T» ptJ U p U a = . (10) 

By computing the corresponding Christoffel symbols for the space-time metric (JT]) and 
substituting their expressions into (fTUj) we obtain that geodesic equations for {/* and 
are satisfied identically, while for the radial geodesic equation we get: 



e 2 * 



0, (11) 



where r is the proper time of the particle and ' denotes partial derivatives with respect to 
r. By considering circular orbits we set again U r = and obtain the following relation: 

- (E - LWf (e- 2 *)' + (0 + 2L (^^) W = 0. (12) 

It turns out that this expression can be easily integrated yielding: 

L 2 {E - LWf 



B 2 e 2 * 



k , (13) 



where k is an arbitrary constant. However, by comparing this result to the relation Q, 
we see that it is in fact a fixed constant: k = — 1. This situation also takes place for the 
null geodesies of photons emitted by stars and detected by an observer located at r — > oo: 
their geodesic equation is equivalent to the normalization condition of their 4-momentum. 

Once we have deleted the kinetic energy component in ([8]) one gets the expression for 
the effective potential on which the star moves ([!]). Thus, by deriving it with respect to the 
radial coordinate and setting the result to zero, we basically impose the minimum condition 
on the effective potential, a necessary condition to have circular orbits for the stars. This 
condition leads to the following expression for E: 

e 23> , L 2 e 2 * 

E = W - V 2L^~' • (14) 
whereas from equation (|9|) we get the following relation: 



E LW +e*^/l + ^ . (15) 
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From both of these equations we see that the second term in the rhs of (|14p must be positive 
definite, leading to the following restriction: 



e 2 * + 



R 2 



^ 0, w's o, 



(16) 



where the last two inequalities must be satisfied simultaneously in the sense that the nu- 
merator and the denominator of that term must have opposite sign for positive L. 

Moreover, from (|14p and (|15p we obtain the expression for the angular momentum L: 



(wy±\w' 













2 -(iy') 2 | 



(17) 



where the following restrictions must be fulfilled in order to have a real angular momentum: 

2R' (e 2 *)' 



i) {W'Y > 



R 3 



-, equivalents, [W') 2 > (e 2 *)' (R)' 



(18) 



this condition makes real the square root of the second term in the radicand's numerator; 

(19) 



ii) (W') ±\W'\\ (W) 2 + 



R 3 

^ \w'\ 



R 2 



^0, 



(20) 



where, again, the last two inequalities must have the same sign simultaneously. 
On the other side, for the energy E we have from (|15|): 



E± 



(W'f±\W'\ \/(W') 2 + ^ w 




21 C 



+ w 



i 



{w ,f ± m ^ {W f ; ^mpL _ <£!^ 









w \ 




2 -(w) 2 | 



(21) 



Moreover, in order for the circular orbits of stars to be stable, we also need the second 
derivative of the effective potential to be positive: 



„2<I> 



1 + 



L 2 
R 2 



+ 2L (E - LW) W" - 2L 2 (W') 2 > 0, 



(22) 



where L and E must be replaced by their expressions (|T7j) and (|2T|) . respectively. 

It is worth noticing that relations (|17p and (|21|) . as well as the condition (|22p . can be 
reduced to the results obtained for the static spherically symmetric approach reported in 
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[TJ [2] when 9 = n/2, W = and R = r under a simple coordinate transformation. For 
instance, by taking N = e* we have: 

2 _ iV 2 2 _ r 3 d r N/N e* rV^ 

£ " i-ra r iv/7v' L " i- r d r N/N and ^"Tt^tF' TT^tF' (23) 

which are precisely the expressions obtained by U. Nucamendi et al. pQ and K. Lake [2], 
respectively, for the energy and angular momentum. 

It is well known that rotation curves of spiral galaxies are inferred from the red and 
blue shifts of the radiation emitted by stars that move in (nearly) circular orbits around the 
central region of the galaxy [3] and that light signals travel on null geodesies with tangent 
4-momentum k^. Here we shall make the assumption that k^ are restricted to lie in the 
galactic plane 9 = 7r/2, and we shall evaluate the frequency shift of a light signal emitted 
from a star in a circular orbit represented by Oe and detected by an observer represented 
by Od- Moreover, we shall suppose that the galactic disk is edge-on directed towards the 
observer, a fact which implies that k° = 00 

The frequency shift associated to the emission and detection of light signals is given by 

l + z = ^ , (24) 

IjJd 

where the frequency of a photon measured by an observer with proper velocity u c \ p c reads 

wc = -k^u£\ Pc , (25) 

and the index C refers to the emission (E) or detection (D) at the corresponding space-time 
point Pc- When the observer is comoving with the particle we have 

ue = ~k^ E . (26) 



Thus, a photon emitted at point Pc in the galactic plane possesses a 4-momentum k c = 
k r , 0, k' p ) c . The corresponding conserved quantities along the null geodesies of light 
signals read 

E 7 = -g^k v = -(g tt k t + g tip k*) , (27) 

£ 7 = g^uVk u = g w k^ + g^k 1 , (28) 
whereas the normalization condition for the 4-momentum k^ka = is 

= gtt {k 1 ) 2 + g„ (k r f + g w (k^) 2 + 2g ttp k t k* . (29) 



and leads to the following expression for the k r in terms of the metric, the conserved 
quantities L 7 , Ely, and k and k v : 

g rr (k r ) 2 = - Lyk v = ^E 2 + 2g tvEl L 1+ g tt Ly ^ 

9 tip ~ 9ttg w 

There are two frequency shifts which correspond to the maximum and minimum values of 

^This fact is also taken into account by astronomers when reporting the corresponding total blue or red shifts 
since they substract the contribution coming from the inclination of the galactic disk from the plane 8 = it/2. 
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uje associated with light propagation in the same and the opposite direction of the motion 
of the signal emitter, in other words, the frequency shifts of a receding and an approaching 
star, respectively. As we shall see later, these max/min values of the frequency shifts are 
reached for stars whose position vector r, with respect to the galactic center, is perpendicular 
to the detector's line of sight, i.e., along the plane where k r = for the observer [2]. 

From the constant character along the geodesies of the product of the Killing vector field 
with a geodesic tangent ([27|) . together with the frequency definition (f25j) and taking into 
account that at (r — > oo), the 4-velocity of the detector is given by = (IP, 0, 0, U v )d, 
we get the following expressions for the detector's frequency 



u D = -fc M u& = - {gttU**? + gt^U'k^ + g^uvtf + g w U^W) \ D = (E^U 1 - L^U*) 
and the emitter's frequency 



D ' 
(31) 



cue = -k^\ E = - (guU'k 1 + g^U**? + g^k 1 + </ w CW + g„U r k r ) \ E 

= (E^U* - L^U* - g rr U r k r )\ E . (32) 

Therefore, with the aid of the frequency shift ()24j) we have for arbitrary star orbits 

_ UE _ {E^-L^-g rr Wk^\ E 
i + Z ~u D ~ {EJJ*-LJJ<P)\ D ■ {66) 

Since the star orbits we are considering are circular, then 

UE _ (E y U t -W)\ E _ Uf s -bU* 

where we have introduced b = -gP- — the impact parameter at infinity, and, hence, |6| 
represents the radial distance at any side of the observed center of the galaxy, where 6 = 0. 
Since we shall consider red/blue shifts either side of the central value b = 0, it is convenient 
to express ([3"4"]) as 

+ £ " {U*-e\b\W)\ D ' (35) 
where e = ±1, and to compute the following quantity: 



rjt ( E-LW \\ 

l + z -ZM.- 1 e2$ )1e (36) 

U D { e 2* )\ D 



where L and E are given by (|17p and (|21[) . respectively. This is the gravitational red shift 
of the center of the galaxy measured by and observer located at r — > oo. Since in order 
to consider red/blue shifts either side of the central value b = 0, we need to subtract this 
quantity from (|35h . we define red and blue shifts as follows 

Z red = z + -z c = v . E . D t rr-ir b , 37 

+ Ujj (Ujj - \b\V%) 11 V ' 

Zblue = Zc-Z-= K f, D f rT-frH , 38 

U f D (Uj) + \b\U%) 11 V ' 

corresponding to receding and approaching stars, respectively. It is obvious that now Z re d / 
Zblue (since z + — z c / z c — Z-). These quantities are the general red and blue shifts, 
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respectively, experimented by light signals traveling along null geodesies and emitted by 
circularly orbiting stars around the center of a galaxy to a distant observer. It is worth 
mentioning that in this formula we have dropped the gravitational red shift of the center of 
the galaxy, a fact that it is indeed taken into account by astronomers when reporting their 
observed data. 

In the special case in which the detector is located "far away" enough from the source of 
information, i.e., if the contribution of the dragging of the detector's inertial frame due to 
the rotation of the system, encoded in Up, is negligible in comparison to the contribution 
coming from the U l D component, in other words, if Up « Up, then the detector can be 
considered static at r — > oo and its 4-velocity will be given by u£ = e 2 *(°°>#. This fact 
can be understood as well from the following analysis: since we are considering the limit in 
which Up << Ujj, and taking into account that Ujj = ^r\ D then 



U%_dip 



= n D « 1, (39) 



where s is the proper time of the orbiting particle (star) and is the angular velocity 
measured by the detector at r — > oo. 

Thus, in the special case in which Up « U l D (or tip « 1,) the general red/blue shifts 
(i37l) and ([38]) reduce to the effective (quasi-static) red/blue shift defined by the following 
expression: 



U 



+ 



71¥- 



Z = "TTtH = (E-LWV J - H ■ ( 4 °) 

U D {.—&*- )\D 

which is symmetric with respect to the center of the galaxy as in [HE]. It is worth noticing 
that even when the dragging of the detector's inertial frame is neglected, the differential 
rotation encoded in W, still plays a nontrivial role in the relation (|39p . Hence we have 



Z 2 



L . (E-LW)W 2 
TPI H Z&6 



Ei3 



I E-LW \ 2 
{ 7** ) 



(41) 



D 



We further need to take into account the light bending due to the gravitational field 
generated by the rotating galaxy, in other words, we need to construct a mapping between 
the impact parameter b and the location of the star r given by its vector position r with 
respect to the center of the galaxy, i.e., the mapping b(r). Following [HE], we shall choose 
the maximum value of Z 2 at a fixed distance from the observed center of the galaxy (at 
a fixed b). From (I4ip it follows that if the function factor that multiplies b 2 is monotone 
decreasing with increasing r, then the maximum observed value of Z 2 corresponds to the 
minimum value of r along the null geodesic of the photons. This minimum value of r 
corresponds to the position of the star either side of the center of the galaxy, lying on the 
plane perpendicular to the detector's line of sight, i.e., on the plane where k r = for the 
observer located at r — > oo. Thus, from (j41|) it follows that the squared impact parameter 
b 2 must also be maximized; this quantity can be calculated from the geodesic equation of 
the photons (or, equivalently, from the k^k^ = relation taking into account that k r = 0) 
and is given by 



-9bp ± J 'dp ~ 9tt9w R 2 W±Re* 9 (R 2 W ± i?e*) 2 

t± = to5 = bl= e« * " 2! 
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Since we look for the maximum value of b, then we should choose either b + or 6_, depending 
on the sign and magnitude of the product RW with respect to e , in such a way that b is 
maximum. Finally, the monotone decreasing condition of the factor that multiplies b 2 in 
(|41|) imposes the following restriction 



, , L (E-LW)W\ 
ln[ R? + ) 



- / 




< 


In ^ 



'E-LW 



(43) 



Thus, the mapping b(r), responsible for the gravitational light bending, is given by (|42p 
under the condition (j4*3j) . 

We should also mention here that the relations (|42p and (|43|) reduce to the results 
obtained by [U [2] when W = and R = r: 



r 



2 



lY / eY 



b =^¥ and (^J < (^^J > ( 44 ) 
respectively. 

As mentioned above, observations are reported as 

Z = v{b) - v{b = 0) 

where Z is given by (|4ip . supplemented by the light bending mapping (|42p . 

If one succeeds in determining the three unknown metric functions <£, i? and W from 
observational data, this would imply that the dynamics of light signals is determined by 
the geodesies of a stationary axisymmetric metric, independently of the assumption of the 
dynamics of the geometry (of a theory of gravitational interactions) or of the nature of dark 
matter, in the case that the latter is needed. However, the task of solving for <£, R and W 
is not trivial at all compared to the spherically symmetric case considered in [TJ [2] where 
there was just one unknown metric function to be determined. 

A way of determining the above mentioned three arbitrary metric functions from obser- 
vations consists in making use of the empiric Persic and Salucci's "universal" formula for 
rotation curves in the halo region of a galaxy [5] : 

" ^ (45) 



1 & 2 + /V 

where the constants a, and correspond to the description of the central, red and blue 
shifts defined in (|36p -( j38j) . and are determined from observations of a given galaxy. An 
alternative way to proceed is to consider just the red and blue shifts, (J3TJ) and (j38|) . related 
to this formula and make use of the relation (|14() (or (I15p if possible) in order to complete the 
system of equations for <£, R and W. Thus, in principle, we can conform a system of three 
nonlinear differential equations of first order for three unknown metric functions. From the 
expressions found for the red/blue shifts it is obvious that the difficulties that arise in solving 
such a system are not trivial at all. Finally, we would like to point out that one can make a 
combined use of the red/blue shifts of galactic rotation curves and gravitational lensing as 
proposed in [191 [20]. However, as stated in [20], till now the combined measured data come 
from different distance scales (red shifts): most high quality rotation curves are available 
for galaxies with a low to intermediate red shift of up to z ~ 0.4, while gravitational lenses 
are easier to detect at intermediate to high red shifts z > 0.4. Therefore, both kinds of data 
are available for the same galaxy, but at different radii and, hence, they are not comparable. 
Thus, it is still difficult to take advantage from both sets of observations simultaneously. 
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This situation will improve in the future when observations with a higher resolution will 
be carried out. In summary, more precise observational data are needed in order to benefit 
from the stationary axisymmetric approach to the galactic rotation curves problem. 
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